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Higher-Order Properties of Analytic Wavelets 
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Abstract — The influence of higher-order wavelet properties 
on the analytic wavelet transform behavior is investigated, and 
wavelet functions offering advantageous performance are identi- 
fied. This is accomplished through detailed investigation of the 
generalized Morse wavelets, a two-parameter family of exactly 
analytic continuous wavelets. The degree of time/frequency local- 
ization, the existence of a mapping between scale and frequency, 
and the bias involved in estimating properties of modulated 
oscillatory signals, are proposed as important considerations. 
Wavelet behavior is found to be strongly impacted by the degree 
of asymmetry of the wavelet in both the frequency and the time 
domain, as quantified by the third central moments. A particular 
subset of the generalized Morse wavelets, recognized as deriving 
from an inhomogeneous Airy function, emerge as having partic- 
ularly desirable properties. These "Airy wavelets" substantially 
outperform the only approximately analytic Morlet wavelets for 
high time localization. Special cases of the generalized Morse 
wavelets are examined, revealing a broad range of behaviors 
which can be matched to the characteristics of a signal. 

Index Terms — Wavelet transform, instantaneous frequency, 
ridge analysis, Hilbert transform, time-frequency analysis 

I. Introduction 

WAVELET analysis is a powerful and popular tool for 
the analysis of nonstationary signals. The wavelet 
transform is a joint function of a time series of interest x(t) and 
an analyzing function or wavelet ij){t). This transform isolates 
signal variability both in time t, and also in "scale" s, by 
rescaling and shifting the analyzing wavelet. The wavelet itself 
can be said to play the role of a lens through which a signal 
is observed, and therefore it is important to understand how 
the wavelet transform depends upon the wavelet properties. 
This permits the identification of wavelets whose higher-order 
properties — with, say, duration held fixed — lead to the most 
accurate representation of the signal. 

Here we focus on analytic, also known as progressive, 
wavelets — complex-valued wavelets with vanishing support on 
the negative frequency axis — defined in continuous time. Such 
wavelets are ideal for the analysis of modulated oscillatory sig- 
nals, since the continuous analytic wavelet transform provides 
an estimate of the instantaneous amplitude and instantaneous 
phase of the signal in the vicinity of each time/scale location 
(t, s). The analytic wavelet transform is the basis for the 
"wavelet ridge" method [1], [2], which recovers time-varying 
estimates of instantaneous amplitude, phase, and frequency 
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of a modulated oscillatory signal from the time/scale plane. 
On the other hand, the analytic wavelet transform can also 
be useful for application to very time-localized structures [3], 
particularly if these features may appear as either locally even 
(symmetric) or locally odd (asymmetric) [4]. The many useful 
features of analytic wavelets are covered in more depth by [5]. 

A promising class of exactly analytic wavelets is the gen- 
eralized Morse wavelet family [6], the joint time/frequency 
localization properties of which were examined by [7]. The 
generalized Morse wavelets have been used to estimate char- 
acteristics of a number of different non-stationary signals, 
including blood-flow data [8], seismic and solar magnetic 
field data [9], neurophysiological time series [10], and free- 
drifting oceanographic float records [11], and have also been 
utilized in image analysis [12]. With two free parameters, 
the generalized Morse wavelets can take on a broad range 
of forms which has not yet been fully explored, and in fact 
this family encompasses most other popular analytic wavelets. 
Two parameters yield a natural characterization of an analytic 
wavelet function, since then the decay at both higher and 
lower frequencies from the center of the pass-band can be 
independently specified. The generalized Morse wavelets will 
therefore be the focus of this study. 

This work has three goals. The first goal, addressed in 
Section 2, is to identify important ways in which higher- 
order properties of analytic wavelets express themselves in 
the wavelet transform. Third-order measures of the degree 
of asymmetry of the wavelet, both in the frequency domain 
and in the time domain, emerge as key quantities relating to 
the precise behavior of the wavelet transform; we note that 
comparable behavior would be expected for discrete analytic 
wavelets at long time scales. The second goal is then to 
establish the particular properties of the generalized Morse 
wavelets, and is accomplished in Section 3. A primary result 
is the existence of a sub-family, characterized by vanishing 
third derivative at the peak of the frequency-domain wavelet, 
which offers attractive behavior for the analysis of oscillatory 
signals. 

The third goal, the focus of Section 4, is to provide practical 
guidelines for the choice of a continuous analytic wavelet 
appropriate for a particular task. When calculating the wavelet 
transform, it is often desirable to choose wavelets to match the 
signal or structure of interest. The generalized Morse wavelets 
emerge as exhibiting a variety of possible behaviors, including 
limits in which the wavelet transform collapses to either the 
analytic filter or, in a certain sense, to the Fourier transform. 
The sub-family mentioned above is shown to derive from 
an inhomogeneous Airy function, and we expect these "Airy 
wavelets" will find value as superior alternatives to the popular 
Morlet wavelet. 

All software related to this paper is distributed as a part of 
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a freely available package of Matlab functions, called JLAB, 
available at the first author's website, http://www.jmlilly.net. 

II. Analytic Wavelet Properties 

In this section, several desirable properties of continuous 
analytic wavelets are introduced: maximization of a conven- 
tional measure of the time/frequency energy concentration; 
the existence of a unique relationship between scale and 
frequency; and minimization of the bias involved in estimating 
properties of oscillatory signals. The second and third of these 
will be shown to relate to the degree of asymmetry of the 
wavelet in the frequency and time domains, respectively, and 
more specifically to third order central moments. 

A. Definitions 

The continuous wavelet transform of a signal x(t) € L 2 (R) 
is a sequence of projections onto rescaled and translated 
versions of an analyzing function or "wavelet" i(>(t), 

f 
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x(u) du 
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where = ip{t)e~ lult dt is the Fourier transform of 

the wavelet and X(u>) is the Fourier transform of the signal. 
Note the choice of 1/s normalization rather than the more 
common 1 / y/s, as we find the former to be more convenient 
for analysis of oscillatory signals. The wavelet is a zero- 
mean function which is assumed to satisfy the admissibility 
condition [13] 
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and is said to be analytic if $>(w) — for oj < 0. The wavelet 
modulus |\I>(w)| obtains a maximum at the peak frequency u^; 
note that w and are both radian frequencies. It will also 
be convenient to adopt the convention that ^(u;,/,)! = 2. 

B. An Overview of Analytic Wavelets 

A commonly used complex-valued wavelet is the Morlet 
wavelet [13], which is essentially a Gaussian envelope modu- 
lated by a complex-valued carrier wave at radian frequency v: 
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As the carrier wave frequency v increases, more oscillations 
fit into the Gaussian window, and the wavelet becomes in- 
creasingly frequency-localized. The second term in (4) and 
(5) is a correction necessary to enforce zero mean, while 
a v normalizes the wavelet amplitude. An expression for a v , 
along with further details of the Morlet wavelet, is given in 
Appendix A. The Morlet wavelet is not, however, exactly 
analytic — it is only approximately analytic for sufficiently 
large v, and this has important implications, as will be shown 
shortly. 



A promising class of analytic wavelets are the generalized 
Morse wavelets, which have the frequency-domain form [7] 



(6) 



where U(ui) is the Heaviside step function and where 

a 0tJ = 2(e 7 //3) /3/7 (7) 

is a normalizing constant. The peak frequency of these 
wavelets is given by = (/3/7) 1 / 7 . The wavelet ipp t7 (t) 
is in fact the lowest-order member of an orthogonal family 
of wavelets for each (/?, 7) pair [7], but we will not be 
concerned with the higher-order members in the present paper. 
The generalized Morse wavelets are the solutions to a joint 
time/frequency localization problem, with analytic expressions 
for both the shape of the concentration region and the frac- 
tional energy concentration [7]. Note that we replace the 
subscript "ip", denoting a property of an arbitrary analytic 
wavelet, with the subscript "/?, 7" for specific properties of 
the generalized Morse wavelets. 

The generalized Morse wavelets form a two-parameter 
family of wavelets, exhibiting an additional degree of freedom 
in comparison with the Morlet wavelet. The nature of this 
additional degree of freedom has not yet been fully explored, 
but it would appear to control variation in higher-order wavelet 
properties with the time and frequency resolution held fixed. 
Because of this adjustability, the generalized Morse wavelets 
can exhibit a very broad variety of behaviors, making (6) a 
fairly general prescription for constructing an exactly analytic 
wavelet. Furthermore, the generalized Morse wavelets sub- 
sume two other classes of commonly-used analytic wavelets, 
the analytic derivative of Gaussian wavelets [3] with 7 = 2, 
and the Cauchy, also known as Klauder, wavelets [13] with 
7=1. 

The number of other existing analytic continuous wavelets 
is fairly limited. There is also the complex Shannon wavelet 
[14, p. 63] 
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where sinc(t) is the sine function and rcct(w) is the unit 
rectangle function. The usefulness of the Shannon wavelet 
ips (t) is limited by its slow (1/t) time decay, a consequence 
of the "sharp edges" of its Fourier transform. Another analytic 
wavelet is the Bessel wavelet [15], defined by 
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where Ki(-) is the first modified Bessel function of the second 
kind. However these two wavelets are less commonly used 
than the Cauchy and Gaussian wavelets that are subsets of 
the generalized Morse family, and in addition, are restricted 
in their behaviors on account of having zero free parameters. 
These considerations justify our focus on the generalized 
Morse wavelets. 
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C. Importance of Analyticity 

The advantage of using precisely, as opposed to approx- 
imately, analytic wavelets such as the generalized Morse 
wavelets was demonstrated by [16], who showed that even 
small amounts of leakage to negative frequencies can result 
in spurious variation of the transform phase. It is important 
to emphasize this point for practical signal analysis. As an 
example, a generalized Morse wavelet and a Morlet wavelet- 
come together with their Wigner-Ville distributions [2], are 
shown in Figure 1. The wavelet Wigner-Ville distribution 
is a fundamental time-frequency object which expresses the 
smoothing implicit in the wavelet transform; see [2] for details. 
The instantaneous frequency [17] — a quantity which reflects 
the time-varying frequency content of a modulated oscillatory 
signal — of each wavelet is also shown; in this case both instan- 
taneous frequencies are very nearly constant, as would be the 
case for a sinusoid. Parameter settings have been chosen such 
that the width of the central time window (as measured by the 
standard deviation of the time-domain wavelet demodulated 
by its peak frequency, defined subsequently) in proportion 
to the period 2it/uj^ is the same for both wavelets. These 
two wavelets appear indistinguishable, and their Wigner-Ville 
distributions are nearly identical. 

However, if we narrow the time window of both wavelets 
by a factor of \/T0, in order to increase time resolution at 
the expense of frequency resolution, we obtain the wavelets 
shown in Figure 2. The Morlet wavelet now exhibits leakage 
to negative frequencies, as well as substantial instantaneous 
frequency fluctuations over the central window. At this very 
narrow parameter setting, a local minimum has developed 
in the amplitude at the wavelet center on account of the 
correction terms. By contrast, the generalized Morse wavelet 
has an instantaneous frequency which is nearly constant over 
the central window, and its Wigner-Ville distribution remains 
entirely concentrated at positive frequencies. Unlike the Morlet 
wavelet, the Morse wavelets remain analytic even for highly 
time-localized parameter settings. This is important for the 
analysis of strongly modulated functions, where the wavelets 
are required to be narrow in time to match the modulation 
timescale. 

The implications of the negative-frequency leakage of the 
Morlet wavelet's Wigner-Ville distribution — a manifestation 
of its departure from analyticity — are drastically degraded 
transform properties. Figure 3 shows the wavelet transform 
of a Gaussian-enveloped chirp using the two wavelets from 
Figure 2. The chirp signal has a frequency which increases at 
a constant rate, passing through zero frequency at time t = 0. 
The negative-frequency leakage of the Morlet wavelet leads 
to interference in the wavelet transform which accounts for its 
irregular structure. Essentially this interference pattern can be 
understood as the interaction of the chirp with its image at 
negative frequencies. Estimates of amplitude or phase proper- 
ties of the signal using the Morlet transform would be badly 
biased. Since the Morse wavelet has no support at negative 
frequencies, the interference is completely suppressed. 

As an aside, we point out that in this example, the derivative 
of the phase of the chirp is shown for reference. However, 
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Fig. 1. A Morlet wavelet (a), a generalized Morse wavelet (b), and their 
respective Wigner-Ville distributions (c,d). These two wavelets are fairly long 
in time, and in a sense that will be made precise later, these two wavelets can 
be said to have the same length. In (a) and (b), the thick solid, thin solid, and 
dashed lines correspond to the magnitude, real part, and imaginary part of the 
time-domain wavelet respectively. In panels (c) and (d), ten logarithmically 
spaced contours are drawn from the maximum value of the distribution to 
1% of that value. The thick dashed lines in (c) and (d) are the wavelet 
instantaneous frequencies. 
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Fig. 2. As with Figure 1, but for parameter settings giving wavelets that are 
very short in time. 



this is not the same as the instantaneous frequency of the 
total signal. The maxima-line of the Morse wavelet transform 
follows the latter rather than the former; see [18] for further 
details. 

This illustrates the importance of analyticity for the de- 
terministic properties of the wavelet transform, the details 
of which have been investigated by other authors. Statistical 
properties are also deteriorated by departures from analyticity: 
whereas an analytic transform applied to a Gaussian process 
produces Gaussian proper transform coefficients [16, p. 424], 
this behavior is lost with a non-analytic analysis function. 
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A Gaussian-Enveloped Chirp 
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Fig. 3. A Gaussian-enveloped chirp is shown in panel (a), while panels (b) 
and (c) show the wavelet transform of the signal with the generalized Morse 
wavelet and Morlet wavelet shown in Figure 2(a) and (b) respectively. In (b) 
and (c), the derivative of the phase of the chirp signal is shown for reference 
as the dashed line. 



Since the advantages of analytic wavelets are well established, 
our concern henceforth will be on the influence of analytic 
wavelet structure on transform properties. We focus on Morse 
wavelets because they are a two parameter family encompass- 
ing most other major analytic wavelets. On account of its non- 
analyticity, the Morlet wavelet is not a valid point comparison 
in most of what follows. We turn now to defining several 
important properties of analytic wavelets. 

D. Mapping Scale to Frequency 

It is common practice to consider the scale s as proportional 
to an inverse frequency. But it is critical to keep in mind that 
any assignment of frequency to scale is an interpretation, and 
there is in fact more than one valid interpretation. The ideal 



wavelet should have these different interpretations be identical, 
such that there is no ambiguity in assigning frequency to a 
given scale. 

One may define three meaningful frequencies associated 
with the wavelet itself: the peak frequency cj^ at which the 
wavelet magnitude is maximized, which is also the 

mode of ^(w)! 2 ; the energy frequency 

'oj\^(uj)\ 2 duj 



oo 
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which is the mean of |\I/(a;)| 2 ; and finally, the time- varying 
instantaneous frequency [17] of the wavelet 



ity(i) = j% {In ^(i)} = ^arg {</>(*)} 



(13) 



evaluated at the wavelet center, u)^,(0). A difference between 
iv^p and uj^ is obviously an expression of frequency-domain 
asymmetry of the wavelet. On the other hand, the energy 
frequency and instantaneous frequency are related by [19] 



(14) 



and therefore a departure of w,/,(0) from uj^ implies that the 
wavelet frequency content is not uniform in time. 

There correspond three separate interpretations of scale as 
frequency. Consider x {t) = cos (uj t), having an analytic 
wavelet transform 



W (t, s) = 1 K) e l 



(15) 



where the contribution from negative frequencies vanishes on 
account of the analyticity of the wavelet. The scale at which 
the magnitude of the wavelet transform is maximum, obtained 
by solving 

J^ (M)| 2 = 0, (16) 

is found to be s = = u^/lo . This is the same as the scale 
at which the rate of change of transform phase is equal to the 
signal frequency, that is, the scale at which 

^{\n[W {t,s)])} = uj (17) 

is satisfied. The peak frequency w$ therefore controls location 
of the amplitude maximum, and the rate of phase progression, 
of an oscillatory feature much broader in time than the wavelet. 
Note that for more general signals (16) and (17) respectively 
define the amplitude and phase ridge curves of the wavelet 
transform [18], from which the instantaneous frequency of the 
signal can be derived. 

Secondly, one may form the energy-mean scale, which 
becomes for the sinusoid 



Si/> 



J™ 8 \W (t, S)\ 2 ds _ J™ 8 (^ Q )| 2 ds 
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(18) 



but following a change of variables, this is seen to be simply 
— uj^/ujo. Thus aty determines the scale at which the 
first moment of the modulus-squared wavelet transform of 
a sinusoidal signal occurs; this gives a integral measure of 
energy content of a signal across all scales. 
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Finally, let Ws (t, s) be the wavelet transform of a Dirac 
delta-function 5(t) located at the origin. The rate of change of 
phase of this transform is 

JU{ln[VMM)])} = ^(£) (19) 

which, at the location of the delta-function, becomes uj^ (0) / s. 
The wavelet central instantaneous frequency lo^ (0) therefore 
controls the rate of phase propagation at the center of a feature 
much narrower than the wavelet. 

We thus have that ui a = ui^/s, uj s = w^/s, and u) s = 
(0) /s define three different mappings of scale to frequency. 
The first will correctly give the frequency of a pure sinusoid 
from the scale s at which its transform obtains a maximum, 
the second will correctly give the frequency of a pure sinusoid 
from the energy-mean scale of the transform, and the third 
fixes the frequency to be the same as the phase progression 
of the transform at the location of an infinitesimally narrow 
impulse. As all of these are mappings arguably correct in 
different senses or for different types of signals, it is desirable 
that all three should be the same. In that case there would be a 
unique and unambiguous interpretation of scale as frequency. 
It will be found that the 7 = 3 generalized Morse wavelets 
have the frequency measures very nearly being equal, while 
maintaining exact analyticity as well as good time localization. 

E. Energy Localization 

The degree of energy localization of a wavelet is conven- 
tionally expressed in terms of its Heisenberg area [2] 



where the standard deviations 

2 ft 2 \m\ 2 dt 
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(22) 



describe the wavelet spread in the time domain and frequency 
domain, respectively; here the standard deviations have been 
defined such that they are nondimensional. The Heisenberg 
area obtains a minimum value of one-half for a function 
which has a Gaussian envelope, but such a function is not a 
wavelet because it is not zero-mean. 

Note that there are other notions of time-frequency local- 
ization. The whole set of generalized Morse wavelets are 
optimally localized in that they maximize the eigenvalues of 
a joint time-frequency localization operator, as shown by [6] 
and investigated in further detail by [7], and indeed this is the 
way the generalized Morse wavelets were initially constructed. 
However, the Heisenberg area is a valuable measure of time- 
frequency localization, since it is in standard usage and permits 
ready comparison among different functions. It will be shown 
later that the 7 = 3 generalized Morse wavelets are close to the 
theoretical minimum of the Heisenberg area, while remaining 
exactly analytic, even for narrow time-domain settings; their 
concentration is comparable to or greater than that of the 
Morlet wavelet. 



F. Minimally Biased Signal Inference 

An important application of analytic wavelet analysis is to 
detect the properties of modulated oscillatory signals of the 
form [1], [2] 



x(t) = a + (t) cos [</>+(£)] • 



(23) 



The amplitude a + (t) > and phase </>+(£) in this model are 
uniquely defined in terms of the analytic signal [20] 
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where U (oj) is again the Heaviside step function. In terms of 
the analytic signal, the original real-valued signal is written as 

x{t) = $t{x+{t)} = 5i{a + (i)e^+ (t) } (25) 

and a + (t) and (f>+(t) are called the canonical amplitude 
and phase [20]. The rates of change of the phase ui(t) = 
d/dt{<p + (t)} and log-amplitude v(t) = d/dt {In [a+(i)]} are 
called the instantaneous frequency and instantaneous band- 
width, respectively. 

Wavelet ridge analysis [1], [2] is a second analysis step 
performed on an analytic wavelet transform which estimates 
the properties of the analytic signal or signals associated with 
the time series. The bias properties of wavelet ridge analysis 
were examined by [18], the results of which we make use of in 
this section. A dimensionless version of the wavelet frequency- 
domain derivative is defined by 

,.n* (n) M 



(26) 



where the superscript "(n)" denotes the nth-order derivative 
of ^>(uj). With this definition, an exact form of the wavelet 
transform of an oscillatory signal, derived in [18], is 



X^(t) EE W (t,LJ^/oj(t)) 

= x+(t)x J 1-1 
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W 3( t ) 
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(27) 



where the quantity 4 (t) is a bounded residual associated 
with truncating the integration outside of a finite range [18]. 

Equation (27) defines a nonlinear smoothing of the analytic 
signal by the wavelet, resulting in a quantity x$ (t) which de- 
parts from the analytic signal x + (i) on account of interactions 
between time-domain derivatives of the signal and frequency- 
domain derivatives of the wavelet. Note that these bias terms 
depend on wavelet derivatives evaluated only at the peak 
frequency ui^. Under a certain smoothness assumption on the 
original signal x(t), expected to hold when the signal is locally 
described as an oscillation, the term in (27) associated with 
is much smaller than that associated with and 
the term associated with ^4(0;^,) — here incorporated into the 
residual — is smaller still; see [18J for details. For this reason it 
is more important to minimize |* 3 ((jj^)| rather than 1^4(0;^ )| 
for a given \f , 2 (w^). 
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Since, as shown below, the square root of normalized 
second derivative 2(^v) * s a nondimensional measure of 
the wavelet duration, it appears desirable choose a wavelet 
with vanishing 4> 3 (lj^) for a fixed value of ^ 2 ( w v)' as tn ^ s 
removes the next-highest-order bias term for a given wavelet 
duration. It will be shown that the generalized Morse wavelets 
achieve this with the choice 7 = 3. We are aware of no 
other analytic wavelets with this property, apart from the 
Shannon wavelets which we do not prefer on account of their 
poor time localization. The Bessel wavelets mentioned above 
have a very large value of the third derivative at the peak 
frequency, ^ 3 (uj^) = —6, which would result in substantial 
bias according to the analysis of [18]. 

III. Properties of Generalized Morse Wavelets 

Having defined three useful properties of analytic 
wavelets — high time/frequency concentration, a unique in- 
terpretation of scale as frequency, and minimized bias for 
analyzing oscillatory signals — our next goal is to form explicit 
expressions of the relevant quantities for the generalized Morse 
wavelets. 

A. Frequency-Domain Moments 

The structure of a wavelet can be described in terms of 
its moments or cumulants. We will use both the frequency- 
domain wavelet moments and the energy moments 

M n ,^ = ±- / w n *(w)dw (28) 

N n ^ = — / w"|*(u;)| 2 dcu (29) 
2tt Jo 

as well as the frequency-domain wavelet cumulants. The 
wavelet moments are the terms in the Taylor series expansion 
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while the coefficients K n; ^ in the expansion 
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(31) 



define the frequency-domain wavelet cumulants. Technically, 
these are called "formal moments" and "formal cumulants" 
since the frequency-domain wavelet is not normalized as a 
probability density function. The wavelet cumulants K n -^ may 
be found in terms of the wavelet moments M n -^ through 
M .^ = exp (i^o;,/,) together with a recursion relation given 
in Appendix B. 

The moments of the generalized Morse wavelets are 

M n -^ = ^(^±11^) (32) 



27i"7 



while the expression 



1 



(33) 



2(l+")/7 " ±n ^^ 

gives the energy moments of the (/?, 7) wavelet in terms of the 
moments of the (2/3, 7) wavelet. Combining (32) with (88) of 



Appendix B, one obtains simple expressions for the cumulants, 
the first three of which are given explicitly in Appendix B. 
In Appendix C it is shown that for the generalized Morse 
wavelets, the series in (30) converges for all t such that \t\ < 1 
for 7 = 1, while for 7 > 1 the radius of convergence is infinite. 

B. Frequency -Domain Derivatives / Time-Domain Moments 

We will also need expressions for the dimensionless 
frequency-domain derivatives ^ n (u) evaluated at the peak 
frequency w^, shown in Section II-F to control the wavelet 
transform of an oscillatory signal. The dimensionless deriva- 
tives can be cast in a form which is somewhat more straightfor- 
ward to interpret. Let the time-domain moments of the wavelet 
demodulated by its peak frequency be denoted by 

/CO 
t n e ~^t^ dt (34) 
-OO 

and recall the correspondence between time-domain moments 
and frequency-domain derivatives [21, Section 5.5], 



*(™)(t, 



(35) 



Since mi ;l /, = 0, the m n -^ are central moments of the demod- 
ulated wavelet e~ luJ ^ t, 4>{t)/mQ- 1 p. This suggests normalizing 
the higher-order demodulate moments by the second moment, 
as one would for a probability density function: 



n;tp 











n/2 ~ % 















n/2 ' 



(36) 



The normalized third and fourth central moments a 3; ^, and 
will be called the demodulate skewness and demodulate 
kurtosis, owing to their formal resemblance to skewness and 
kurtosis of a probability density function. However, it is 
important to keep in mind that the demodulated wavelet is not 
a probability density function since it is in general complex- 
valued and not nonnegative. 

Now define a dimensionless measure of the wavelet time- 
domain length or duration 



2 / m 2;j, 

V" — 77 ~o — / = ^ 



27:/^ 



717-2;^) 
»TlO ; V> 



(37) 



such that P^/tt is the number of oscillations at the peak 
frequency which fit within the central wavelet window, as mea- 
sured by the standard deviation of the demodulated wavelet. 
Increasing increases the frequency-domain curvature in the 
vicinity of the peak frequency, narrowing the wavelet in the 
frequency domain and hence broadening the wavelet in the 
time domain. 

For the generalized Morse wavelets, we find in Appendix D 
that P/3 j7 is simply y/Jh/, while the demodulate skewness and 
kurtosis become 

.7 — 3 .7 — 3 



0,1 



2 2 

= 3-3{a 3;/ 3, 7 } --753- 



(38) 



(39) 
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where we note the demodulate skewness is a purely imaginary 
quantity. Note that thej;hoice 7 = 3 causes the demodulate 
skewness [and hence ^ 3 (w^)] to vanish — while simultane- 
ously maximizing the magnitude of the demodulate kurtosis 
«4 ; /3,7 for a fixed value of Pp a . As mentioned in Section II-F, 
we are less concerned with the kurtosis than with the skew- 
ness: for signals which are locally oscillatory, it is expected 
that nonzero skewness will contribute more substantially to a 
bias of (27), the analytic signal estimated from the wavelet 
transform [18]. 

The dimensionless duration Pp a and demodulate skewness 
^{®3;p,~/} form a natural two-parameter description of the 
generalized Morse wavelets. Particular values of Pp a and 
^{ a 3-,p,-y} give a unique (/?, 7) pair. Pp tl is a normalized 
second-order moment, while a^p n is a normalized third- 
order moment measuring the degree of asymmetry of the 
demodulated wavelet in the time domain. For a given duration 
Pp,-y, a range of shapes can be obtained for the generalized 
Morse wavelets by adjusting 3{a3 ;( 3, 7 }. The fourth-order 
behavior, as expressed by the demodulate kurtosis a4-,p,j, is 
not free, but is an implicit function of the two lower-order 
quantities Pp n and 9{a3 ;/ 3 ;7 }. This clarifies how (3 and 7 
translate directly into controlling the wavelet moments. 

The generalized Morse wavelet parameters 7 and [3 are 
plotted as a function of Pp n and 3 {&3;p, 7 } in Figure 4(a), (b). 
The parameter j3 increases with increasing Pp tl but decreases 
with increasing 3{a3 ; /3 i7 }. The 7 curves, on the other hand, 
change character at 7 = 3, with contours of lower values being 
concave down and those of higher values being concave up; 
7 = 3 itself is the horizontal line 3 {a^-p^} = 0. For a given 
Pp tl with j3 > and 7 > 0, ^{a^p^} is bounded from 
below by —3/Pp n for 7 = 0, but has no upper bound; this 
lower bound is the cause of the empty region below the 7 = 
contour in these plots. 

C. Wavelet Frequency Measures 

We now return to the question of assigning a frequency 
interpretation to the wavelets. From the results of the preceding 
sections, we see that the energy frequency defined in (12) 
becomes 

while the time-varying instantaneous wavelet frequency (19) 
is [using (31)] 

<M*) = j t ^{\^{t)}=K 1 .^-^K^t 2 + ...m 

when expanded in terms of the wavelet frequency-domain 
cumulants; only the odd cumulants appear because of taking 
the imaginary part. Thus = K\.$ = Mi ; ^/M 0; ^. The 

quantity 



1 d 2 uj^ _ K 3 .^ 

2;V 



2;i/> 



K. 



(42) 



is a nondimensional measure of the curvature of the in- 
stantaneous frequency evaluated at the wavelet center. The 



instantaneous frequency curvature has a simple interpreta- 
tion for a real-valued, nonnegative definite frequency domain 
wavelet such as a generalized Morse wavelet. For such a 
wavelet, $(w)/Mo ; ^, is a probability density function having 
K^ip/ K^. 2 as its coefficient of skewness, which is the neg- 
ative of the dimensionless instantaneous frequency curvature. 
Thus frequency-domain skewness corresponds to curvature of 
the instantaneous frequency. 

One finds for the generalized Morse wavelets that the energy 
frequency is 



N 



1 M 1;2 p, 



p { 2/3+2 
1 I 7 



iVo ;/ 3, 7 2V7 M 0;2/ 3, 7 2V7 p 



(43) 



while 



W/3, 7 (0) = Kx-p, 



M x .p, 
M -p. 



r(^) 



2 l/ 7 



(44) 

gives the time-varying instantaneous frequency at the wavelet 
center. An expression for the frequency curvature may be 
found from (32) together with (85) and (87) of Appendix B. 

The behaviors of these frequency measures of the general- 
ized Morse wavelets versus Pp^/ir and 3 {a 3 -p n } are shown 
in Figure 4(d),(e),(f). A change in character is observed at 
7 = 3. At 7 = 3, the ratio of the energy frequency to the 
peak frequency u)p^/u)p t7 , and the ratio of the instantaneous 
frequency at the wavelet center (t = 0) to the peak frequency 
u)p^(Q)/cjp^, are both very close to unity, except for very 
short wavelet durations with Pp^/n < 1. For larger values of 
7, or positive 9 {a 3 -p^}, both of these two ratios are generally 
smaller than unity, while for smaller values of 7 or negative 
^{ a 3-,p,-f} both ratios generally exceed unity. The exception 
is for very short durations Pp, 7 /n < 1, where one finds these 
ratios becoming increasingly large and positive as Pp^jit 
decreases with fixed 3 {a3 ;( g, 7 }. Meanwhile the instantaneous 
frequency curvature at the wavelet center, Figure 4(f), exhibits 
a similar pattern but with the sign reversed. The 7 = 3 
wavelets have negligible instantaneous frequency curvature 
except as the duration becomes very short. 

For Pp n /n > 1, the 7 = 3 wavelets have ujp n w wp n w 
ujp n (0) to a very good approximation, nearly obtaining the 
unambiguous interpretation of scale as frequency which was 
sought in Section II-D while remaining exactly analytic. That 
utp.y and cop.-y should be almost identical for some value of 
7 is not obvious. The former is a simple algebraic expression 
in terms of powers of (i and 7, whereas the latter is given by 
a ratio of gamma functions. In Appendix E it is shown that 
for increasing /?, the ratio of these two quantities converges 
rapidly to unity for 7 = 3 due to the asymptotic behavior of 
the gamma function. 

The sign change of the wavelet frequency curvature ob- 
served in Figure 4(f) gives the border between two qual- 
itatively different behaviors, as is seen in Figure 5. For 
negative curvature, we have concave wavelets in which the 
instantaneous frequency takes on its maximum value at the 
wavelet center, while for positive curvature the wavelets are 
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Properties of Generalized Morse Wavelets 




0.5 1 2 4 8 16 0.5 1 2 4 8 16 0.5 1 2 4 8 16 

Duration P„ In Duration P D In Duration P. In 

Por P.y P.y 



Fig. 4. The behavior of Morse wavelet quantities as a function of Pp n and the imaginary part of demodulate skewness a3 ;/ 3, 7 is presented. All panels have 
the same axes, with Pp^/n being the x-axis and 3?{a 3;( 3 7 } being the y-axis. The six panels show six different quantities contoured as a function of this 
plane. Panels (a) and (b) show 7 and /3, respectively. Heavy solid lines show 7 = 1 and /3 = 1, white lines with black outlines show 7 and /3 = 2, 3, and 4, 
and thin solid lines show 7 and f3 = n 2 for integer n with 3 < n < 10. The dashed line in all panels is the 7 = contour. The Heisenberg area is shown 
in panel (c) with a contour interval of 0.01 from 0.51 to 0.59; the heavy solid line is the 0.51 contour. Panel (d) shows uip^/uip^ — 1, the difference of the 
ratio of the energy frequency to the peak frequency from unity, and similarly panel (e) shows u/g t7 (0)/uig t ^ — 1 where ii^ 7 (0) is the value of the wavelet 
instantaneous frequency at the wavelet center. Finally panel (f) gives the dimensionless curvature of instantaneous frequency ibp^(t) as defined in (42). The 
last three panels all have the same contours, which range from -0.2 to 0.2 with a contour integral of 0.025. Thin solid contours are for positive values, white 
contours with black outlines are for negative values, and the heavy curve is for the zero contour. Note that 7 = 3 lies along the x-axis. 



convex and the instantaneous frequency takes on a minimum 
value at the wavelet center. Regions of large amplitude thus 
correspond to regions of high frequency for the concave case, 
but to regions of low frequency for the convex case. The 
degree of convexity, or concavity, controls how the wavelet 
filter will respond preferentially to signals having frequency 
minima, or maxima, at the wavelet center. For Pp.-y/n > 1, the 
7 = 3 wavelets are the division between these two cases. The 
7 = 3 wavelets have very small instantaneous frequency cur- 
vature, and have Wigner-Ville distributions which are roughly 
symmetric about the central frequency, as in Figure 5(b). It 
was seen earlier in Figure 2(d) that this symmetry becomes 
compromised for very time-localized settings, corresponding 
to "squashed" appearance of the Wigner-Ville distribution and 
to the curvature apparent in Figure 4(f) for small Pp^/n. 



D. Energy Localization 

We now address the problem of time/frequency concentra- 
tion as measured by the Heisenberg area. The frequency spread 



defined in (22) simplifies to 



2 



N- 



2;V ~2 



N 0;i> 



(45) 



making use of the definition of uij, (40). For the time-domain 
spread (21), note that one may show 



t;V> 



= UJ 



(46) 



using the relation between time-moments and frequency- 
domain derivatives together with Parseval's theorem. This 
becomes for the generalized Morse wavelets 



; *;/3,7 



73,7 



•Wu;/3,7 



, N, 



0;/3-l,7 



a 



+ 7 



2-^0;/3-l+7,7 



/3-1.7 



a 



N, 



0-1+7,7 



O;0-l+7/2,7 
,2 

^-1+7/2,7 



(47) 



which can then be expressed using (7) and (33). 

The Heisenberg area is shown in Figure 4(c) as a function 
of the wavelet duration Pp :7 /n and demodulate skewness 
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S { a 3-,p,i}- It is clear the wavelets exhibiting small time- 
domain skewness have a small Heisenberg area. As Pp, 7 /iT 
increases, the Heisenberg area approaches the limiting value of 
one-half for any value of the demodulate skewness 3 {a^-^^}. 
The theoretical minimum value of the Heisenberg area of 
one-half is not quite obtained, except in the limit of long 
duration, evidently on account of asymmetry induced by the 
constraint of analyticity. However, we may point out that 
perfect concentration is not achieved by the Morlet wavelet 
either — although constructed from a Gaussian, the existence 
of the correction terms lead to departure from the theoretical 
minimum value of one-half. In fact, numerical computations 
we have performed (not presented here) show the Heisenberg 
area of the 7 = 3 wavelet is comparable to or smaller than 
that of the Morlet wavelet. 

Summarizing the results in this section, we see from Fig- 
ure 4 the special properties of the 7 = 3 wavelets. In the 
vicinity of 7 = 3, the generalized Morse wavelets obtain their 
minimum Heisenberg area, and for sufficiently large duration 
Pp n the peak frequency W/3, 7 , energy frequency and 
central instantaneous frequency &p n {Q) all become indistin- 
guishable while the wavelet instantaneous frequency curvature 
(42) vanishes. The 7 = 3 generalized Morse wavelets therefore 
obtain the ideal behavior with respect to the three issues 
raised in the Section 2 for Pp^/n > 1. For extremely 
time-concentrated wavelets with Pp^/ir < 1, the curves 
along which the frequency pairs are identical, the frequency 
curvature vanishes, and the Heisenberg area is minimized, all 
begin to diverge from one another. 

IV. Special Cases of Generalized Morse Wavelets 

In this section we step back from an emphasis on prop- 
erties relevant for analysis of oscillatory signals. Instead we 
examine the broad variety of behavior of the generalized 
Morse wavelets, with the idea in mind that these could be 
considered a generic family of analytic wavelets appropriate 
for analyzing many different types of signals. In this section 
we therefore isolate special cases of these wavelets, exploring 
the boundaries of the family as well as the relationships among 
its different members. 



A. Interpretations of (3 and 7 

We discuss two important interpretations of (3 and 7, the 
first pertaining to the relationships among different members 
of the generalized Morse wavelet family, and the second to the 
time-domain and frequency-domain decay of the wavelets. 

1) Differentiation and Warping: The (/3, 7) generalized 
Morse wavelet was represented as a nonlinear transformation 
of the {^f- — 5,1) wavelet by [7]. This representation 
was crucial for deriving the localization properties of the 
generalized Morse wavelets in the time-frequency plane for 



(3 > > 0, since the time-frequency localization operator 
for which these wavelets form the eigenvectors [6] is only well 
defined in this case. Here we present an alternate construction 
which more directly reflects the different roles of the (3 and 7 
parameters. 



Note that over the entire range 7 > and (3 > 0, the inverse 
Fourier transform ip/3 n {t) of (6) still defines a valid filtering 
function. We refer to ipp n (t) over this entire range as the 
generalized Morse filter, only a subset of which corresponds 
to the generalized Morse wavelets. In fact, it is easy to see 
that this filter is zero-mean for (3 > 0. Computing defined 
in (3) for the generalized Morse wavelets, we find [using (32)] 



C/3, 7 



7^22/3/7+1 



7 



(48) 



and over the entire range 7 > and (3 > admissibility is 
satisfied, as is the constraint of finite energy. Thus only (3 = 
or 7 = are not valid wavelets. 

Note that the time-domain (3 = 0, 7 > 1 filter 



^0,7 (*) 



1 r 

n J-< 



(49) 



may be expressed in terms of the frequency-domain (3 = 0, 
7=1 filter as 



^0 



, 7 (*) = ^ f^ ' 1 ^' 



(50) 



since ^ ,i(w) = 2e~". Equation (50) states that the 
frequency-domain power distribution of ^0.1 ( w ) is mapped 
onto different Fourier components through the substitution 
to 1 > lo^ . But substituting the inverse Fourier transform, 

\&o,i(w) = ipo,i(t)e~ luJt dt, one may write instead 



^0,7 (*) 



where we have defined 



f 

j — i 



ipo i(u)K^(t, u) du 



■dw 



2?r Jo 



(51) 



(52) 



as a time-domain transformation kernel function. Note that 
for 7 = 1 one has K\{t,u) = S(t — u) where 5(t) is the 
Dirac delta-function. Thus incrementing 7 is accomplished by 
a frequency-domain warping. 

Subsequently, the (3 > 1 filter is obtained from the (3 = 
filter for fixed 7 and (3 € N via the time-domain differentiation 



#,7(0 



, Nfl 1 d p , , . 



(53) 



Therefore all generalized Morse wavelets can be generated by 
beginning with V'o.i making this function analytic to 
obtain ipo.i(t), warping the frequency content to increment 7, 
and then differentiating in the time domain to increment (3. 
The nature of the originating function 5i{V'o,i(t)} will De seen 
shortly. 

2) Frequency and Time Decay: Clearly the parameter 7 
also controls the high-frequency decay of the wavelet. We now 
show that (3 controls the time-domain decay. The time-domain 
form of the generalized Morse wavelets is expressed by the 
inverse Fourier transform 



#,7(0 



1 r 
2^ Jo 



duj. 



(54) 
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Morse with y=2 and (3=6 Morse with 7=3 and |3=4 Morse with y=4 and (3=3 




-150-100 -50 50 100 150 -150-100 -50 50 100 150 -150-100 -50 50 100 150 
Time Time Time 

Fig. 5. The time-domain generalized Morse wavelets (a,b,c) and their Wigner-Ville distributions (d,e,f). All three wavelets have the same value of 
Pg^{ujp tl ) = y/j3j = 2\/3, but differing values of 7 as indicated in the captions. Line styles and contour intervals are as in Figure 2. 



One may obtain their asymptotic time domain behavior using 
the method of [22, p. 407] by noting that 



ui' e 



s=0 



(55) 



Inserting this into (54), we find that the integrals of the terms in 
this summation, while possibly divergent, are Abel summable 
[22, p. 407] and it follows from this reference that 



V>/3, 7 (t) = a/?, 



E 

s=0 



H)!exp|^ + /3+1 >>x 



r(s 7 + /? + 1) 



We therefore obtain the asymptotic behavior 



^18,7 (*) 



m(/3+n/2 r (ff+ 1 ) 
«/3,7 e t p+l : 



t 



OO 



(56) 



(57) 



since the smallest power of 1/t dominates at large times. The 
0{t~^ +1 ^) behavior could have been anticipated from the 
fact that the frequency-domain wavelet ^^^(u) is (3 times 
differentiable but has a singularity in the ((3 + l)st derivative 
at lo = 0. 



The different roles of (3 and 7 are illustrated in Figure 6, 
which shows the first sixteen generalized Morse wavelet filters 
at integer (3 > and 7 > 1. Note that the appearance of 
the twelve wavelets ((3 7^ 0) varies dramatically despite the 
fact that all have been set to have the same peak frequency. 
The action of differentiation (increasing (3) is to broaden the 
central portion of the filter, while at the same time making the 
long-time decay more rapid. On the other hand, increasing 
7 reduces the curvature of the filter envelope at its center, 
also causing it to broaden, but without changing the long-time 
decay. This broadening of the central window width as (3 or 7 
increases agrees with our earlier identification of P ( g 7 = y/JJj 
as a dimensionless time-domain duration. Adjusting (3 and 7 
together therefore permits the "inner" width of the wavelet 
window to be controlled independently from the long-time 
decay. 

3) Symmetry Versus Compactness: Earlier it was shown 
that time-domain symmetry of the demodulated wavelet is 
controlled by 7 through the "demodulate skewness" parameter 
a 3;/3,7 — i(l ~ 3)/Pg 7 (38). Thus we can interpret (3 as 
the decay, or compactness, parameter and 7 as the symmetry 
parameter. Note the differing behaviors of the wavelet with 
fixed Pfj tl = y/Jhj on either side of 7 = 3. For 7 < 3, 
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Cauchy Family (7=1) Gaussian Family (y=2) Airy Family (7=3) Hyper-Gaussian Family (7=4) Frequency Domain 




Fig. 6. The time-domain forms of the generalized Morse filters for 7 = 1^1 and /3 = 0-3 are shown in the first four columns, while the fifth column shows 
the frequency-domain version of filters with 7 = 1-4 for each value of (i, with line styles as labeled. All wavelets with /3 7^ have the time axis normalized 
by Pp n , while the (3 = filters have been rescaled along the time axis such that their first time-domain derivatives (multiplied by a constant) are the (3 = 1 
wavelets shown. 



time decay increases as (3 increases from a minimum at 
(3 = Pp 7 /3, and the corresponding decrease in 7 makes 
the wavelet less symmetric. On the other hand, for 7 > 3, 
decreasing (3 from a maximum at /3 = Pj /3 also makes the 
demodulated wavelet less symmetric. In this case the wavelet 
is most symmetric when its time decay is also strongest, and 
this occurs at 7 = 3. Time-domain symmetry and compactness 
are therefore antagonistic for 7 < 3 but covary for 7 > 3. 

B. Domain Boundaries 

Next we examine the behaviors of the generalized Morse 
wavelets at extreme values on the (/?, 7) plane. 

1 ) The Analytic Filter Family: The frequency-domain gen- 
eralized Morse filter for (3 = 7 = is simply twice the unit 
step function 

*o,oM = 2U(u) (58) 
so that the time-domain Morse filter is the analytic filter [23] 

Mt) = Ht) + —+ (59) 

which, of course, is not a wavelet. Application of this filter 
to a signal x{t) simply recovers the analytic version of the 
signal, i.e. 

W x , 0fi (t,s) = x+(t), (60) 



independent of scale. This shows that the generalized Morse 
filter includes the analytic filter as a special case. Commuting 
the differentiation operator with the analytic filter in (53) 
shows that taking the wavelet transform with the wavelet 
^3,0 (*) for /3 G Z essentially involves taking the /3th derivative 
of the analytic signal. 

2) The Complex Exponential Limit: The generalized Morse 
wavelets have an interesting behavior in the case (3 — ► 00. It 
follows from asymptotic expression for gamma function ratio 
(102) in Appendix E that 

MpM „ WW* = /?^oo. (61) 

Now, the moment expansion of the wavelet (30) can be 
rewritten as 

* 9 „( t >/Mo,,, - E<f (62) 

n=0 u ' p ' 7 

while a complex sinusoid at the wavelet peak frequency has a 
Taylor-series expansion 

e ^, Tt = J-^-lco^r. (63) 

n=0 

From (61) we then see that for any fixed n, the nth moment 
of the normalized wavelet ifjp n (t)/M -/3 :7 becomes identical 
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with the nth moment of the complex sinusoid e 1 



as using (53) and where 



(3 approaches infinity. In this sense the generalized Morse 
wavelets approach a complex sinusoid as (3 increases. Equating 
(62) and (63) over some range of times would however 
require careful consideration of terms with n = 0(f3) in the 
summations. 

C. The 7 Families 

Finally we examine in more detail the generalized Morse 
wavelet families for the first few integer values of 7. 

1 ) The Cauchy Wavelets : The 7 = 1 family corresponds to 
the Cauchy wavelets [13, p. 28-29]. For (3 > and 7 = 1 the 
generalized Morse filter becomes the analytic Cauchy filter 

1 f°° 1 



+ i- 



dui 
i 



7r(l - it) 



(64) 



(65) 



TT(l+t 2 ) 7T(l+i 2 ) 

such that 3? {-00,1 (01 * s tne Witch of Agnesi curve, or, 
equivalently, the standard Cauchy probability distribution. This 
filter therefore specifies the joint effect of applying the analytic 
filter and smoothing by the Witch of Agnesi. The /3th filter for 
f3 e N and j3 > 1 is then obtained by (53) to be 



(e/Pf - 



1 r(/3 + i) 



(66) 



7T (1 - ' 

a form which, it turns out, is in fact valid for all (3 > [7], 
not just the integers. 

Following the results of Section IV-A1, all generalized 
Morse filters with (3 e N and 7 > 1 are generated from the 
Witch of Agnesi (the real-valued curve in the upper of left- 
hand corner of Figure 6) by analytization followed by warping 
followed by differentiation. Although not itself a wavelet, this 
time- and frequency-localized function forms the basis for all 
generalized Morse wavelets, and can therefore be thought of as 
the "queen mother wavelet" function. The next two subsections 
demonstrate that the 7 = 2 and 7 = 3 warpings generate 
two other important functions, the Gaussian probability density 
function and the inhomogeneous Airy function. 

2) The Analytic Derivative of Gaussian Wavelets : The 
7 = 2 family corresponds to analytic Derivative of Gaussian 
wavelets [3]. With (3 = and 7 = 2, the generalized Morse 
filter becomes 



^0,2 (£) 



= l f 

n Jo 



cLu 



2^ 



~t 2 /4 



+ i- 



D(t/2) 



(67) 



(68) 



where D (t) = e~ l J a e u du is the Dawson function. This 
extends the representation of [7, p. 2667] which is only valid 
for the real part and for even values of j3. The ((3, 2) wavelets 
are given for integer (3 > 0, with Hp{x) denoting the /3th 
Hermite polynomial [24, eqn. 22.2.14], by 



\2 



H P (t/2) e- t2 / 4 + i(-lf4= D {f3 \t/2) 



(69) 



D {n) {t) = (-l)"x 

J H n (t) D(t) -J2(k) Hn - k {t) lk ~ lRk -^ W \ (70) 



k fe=i 
gives the form of the nth derivative of the Dawson func- 
tion, which has been derived using Leibniz's theorem 
[24, eqn. 3.3.8]. 

These wavelets have been proposed for singularity analysis 
by [3], but no analytic expression for their time-domain form 
has been given previously as far as the authors are aware. As 
illustrated in Figure 5(d), the instantaneous frequency curve 
for the analytic Derivative of Gaussian wavelets is concave; 
this is true for all (3 as Figure 4(c) shows. Their frequency 
domain behavior makes them less appropriate for the analysis 
of oscillations than the 7 = 3 wavelets. 

3) The Airy Wavelets : The 7 = 3 generalized Morse 
wavelets in fact derive from an inhomogeneous Airy function, 
therefore we suggest calling this family the Airy wavelets. 
The second inhomogeneous Airy function Hi(z), also known 
as the second Scorer function, is defined by the integral 
[24, p. 448, eqn. 10.4.44] 



m(z) 



1 r 

n Jo 



e- u3/3 e zu du. 



(71) 



Thus the generalized Morse filter with (3 = and 7 = 3 is 
simply 



^0,3 (*) - 



1 



dw = 



3V3 



Hi 



it 
3V3 



(72) 



which is the inhomogeneous Airy function evaluated at an 
imaginary argument. Differentiating the analytic Airy filter 
00,3 (t) (3 times, as in (53), one obtains 

, , s , -x/3 1 1 d 



2 3V3 



Hi 



it 
3i73 



(73) 



as an expression for the /3th Airy wavelet ipp^if) with (3 E N 
and (3 > 1. Note that the (3=1 Airy wavelet is not within 
the localization regime (3 > (7 — l)/2. Examples are shown 
in Figure 2(b) and Figure 5(b). As already discussed, the 
instantaneous frequency of the wavelet is nearly constant over 
the width of the wavelet, and the wavelet function exhibits no 
preference for its Wigner-Ville distribution to skew to smaller 
or larger frequencies on its periphery. 

4) The Hypergaussian Wavelets : The 7 = 4 family does 
not have an analytic time-domain expression in terms of 
known functions as far as the authors are aware. However, this 
family is interesting because it is the first integer 7 family ex- 
hibiting convex behavior of the instantaneous frequency curve. 
We may note that the analytic Gaussian filter is generated from 
the analytic Cauchy filter via the frequency-domain warping 

Mt) = tt / *o,i(wV wt du; (74) 

271" J-00 

while the analytic 7 = 4 filter may be expressed as 

iM*) - 7T / *o,2(w 2 )e^du>. (75) 
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Thus the relation between tpo^t) and '00,2 (*) is the same 
as that between Vo,2(i) and tpo,i{t)- We therefore suggest 
"Hypergaussian" as a name for ^0,4 (*) since it involves a 
second iteration of the nonlinear operation creating the analytic 
Gaussian filter from the analytic Cauchy filter. 

V. Discussion and Conclusions 

This paper has examined the higher-order properties of 
analytic wavelets and their impact on the behavior of the 
wavelet transform. Three important wavelet properties were 
discussed — time-frequency localization in terms of the Heisen- 
berg area, the existence of a unique correspondence between 
scale and frequency, and minimized bias in the extraction of 
oscillatory signals. The latter two were shown to be related 
to third order moments of the wavelet. These properties were 
examined for the generalized Morse wavelets, a two-parameter 
family of exactly analytic wavelets. 

The existence of a unique correspondence between scale 
and frequency requires symmetry about the peak frequency, 
as measured by the frequency-domain skewness, and also 
equality between the mean and the mode of the squared 
modulus of the frequency-domain wavelet. Minimized bias in 
estimating instantaneous properties of modulated oscillatory 
signals was found to require a vanishing third derivative at 
the wavelet peak frequency, which is equivalent to a vanishing 
third central moment of the time-domain demodulated wavelet. 
Thus with a lower-order property held fixed — such as the 
wavelet duration in proportion to its period, denoted here by 
P$ — choosing a wavelet which has a high degree of symmetry 
in both the time domain and the frequency domain leads to 
good performance for the analysis of oscillatory signals. These 
results for continuous analytic wavelets could also contribute 
to an improved understanding of the behavior of discrete 
analytic wavelets, such as those of [5]. 

One member of the generalized Morse wavelet family was 
found to have zero asymmetry in the time domain, as measured 
by the third central moment of a demodulated version of itself, 
as well as competitive performance in terms of other criteria. 
This is the 7 = 3 wavelet, shown herein to be derived from 
an inhomogeneous Airy function. In fact the Airy wavelet 
preserves the spirit of the Morlet wavelet more than the Morlet 
itself, remaining nearly symmetric in the frequency domain 
and maintaining a nearly optimal Heisenberg area even at 
high time concentration, yet without compromising its exact 
analyticity. 

The roles of the two parameters 7 and j3 in setting practical 
properties of the wavelet filters was investigated in detail. 
Here we showed that the former controls the width of the 
inner wavelet window without impacting the time decay, while 
increasing the latter broadens the wavelet central window but 
increases the rate of decay at large times. The generalized 
Morse wavelets include as special cases the Cauchy wavelets 
(7 = 1) as well as analytic versions of the Derivative of 
Gaussian wavelets (7 = 2). The Airy wavelets emerge as 
the approximate boundary between two qualitatively different 
sorts of behavior, which we identify as convex or concave 
depending upon the sense of curvature of instantaneous fre- 
quency curve. The broad range of behavior of the generalized 



Morse wavelets, together with their attractive properties for 
certain values of j3 and 7, suggests their use as a generic 
family of exactly analytic wavelets. 

Appendix A 
The Morlet Wavelet 

In this section we address some details of the Morlet wavelet 
(4-5). First we find an expression for the peak frequency uj v 
at which the frequency-domain wavelet obtains its maximum 
value, which is not the same as the carrier frequency v. The 
peak frequency w„ of the Morlet wavelet occurs where the 
first derivative 

*'„(w) = a v e-* { - u - v ?[v + u){e- uv -\)] (76) 

vanishes. This occurs when 

w-v = ue^ u (77) 

the solution to which may be found by introducing v = v/u v . 
Then (77) leads to 



u v {v) = \ ~ 



In (1 - v) 



(78) 



and since < v < u v , we can numerically solve (78) on the 
interval < v < 1 to obtain u v (v). Setting 



v 



(79) 



for the normalization function a v obtains our chosen value of 
$„(w„) = 2, and from the above parametric form for lo v we 
likewise know a v as a function of the carrier wave frequency. 

Additional wavelet properties are given by the value of 
higher-order derivatives at the peak frequency. One may verify 

= - [w(w - v) + 1] tf„(w) - [2w - v] tf'„(w) (80) 

as an expression for the second derivative in the frequency 
domain, which leads to 



*< 2) (w„) = -ul[w v (u v -v) + l] 



(81) 



for the normalized second derivative evaluated at the peak 
frequency. The wavelet duration is then 



P v = \J -^\bJ v ) = W v y/u v {u}„ -v) + l 



(82) 



and one may note that as v becomes large, one has ld v <~ v, 

^>i?\uj„) <~ -v 2 , and P v ~ v. 

Appendix B 
Wavelet Moments and Cumulants 

To find the relation between the wavelet cumulants and the 
moments, note that 



i>{t) = exp (In [iP(t)}) = exp ^ ^-K^ 

\n=o n - J 



71=1 



n) 



-B n {Ki-rf, K 2 -^, ■ ■ ■ K n .^) 



(83) 



[using (31)] which implicitly defines B n (ci, C2, . . . c„), the 
nth-order complete Bell polynomial; see [18] and references 
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therein for details. Then equating powers of t between (30) 
and (83), one finds the moments are given in terms of the 
cumulants as M 0; ^ = cxp (ATo;v) an d 

^ = B n {K h ^K 2 ,^...K n ,^) n>\. (84) 

Inverting (84) leads to 

M i;/3,7 



Ki-^j, — 



K 



2;l/> — 



M 0;/3, 7 



M 0; 



M, 



(85) 
(86) 



0;j/> 



as expressions for the first three cumulants. These differ from 
the usual expressions between moments and cumulants [e.g. 
25] because the frequency-domain wavelet is not normalized 
as a probability density function, that is, Mq-^, ^ 1. More 
generally, the recursion relation 



n;tp 



£7 I 7v, 



fc=l 



k - 1 



Mi 



0;if> 



relates the moments to the cumulants and vice-versa. 

Appendix C 
Convergence of Morse Moment Expansion 

Here we investigate the convergence of the moment ex- 
pansion (30) for the generalized Morse wavelets. The series 
converges for all t such that \t\ < r, where r is a positive 
constant referred to as the radius of convergence [26, p. 203]. 
This radius may be determined by the ratio test as 

nl M n+1 .^ 



= lim 



n — >oo (n + 1)! M„ ;1 /, 



(89) 



and one finds for the generalized Morse wavelets with fixed 
09,7) that 



= lim 



1 



+oo n + 1 p 

lim n 1 ^- 1 {l/ 1 ) lh 

I >OQ 



(90) 



(91) 



using (32) and the asymptotic behavior of the gamma function 
given subsequently in (102). The moment expansion for the 
generalized Morse wavelets therefore has radius of conver- 
gence r = 1 for 7 = 1, and infinite radius of convergence for 
7 > 1. 

Appendix D 
Wavelet Frequency-Domain Derivatives 

To find the generalized Morse wavelet frequency-domain 
derivatives, first note that there exists a simple expression 



ln*^) = [(-lr-^n - 1)!] - W J] (7 - p) 

p=0 

(92) 



[n > 1] for the derivative of the logarithm of the wavelet. 
Taylor-expanding the frequency-domain generalized Morse 
wavelet about any fixed frequency uo a leads to 



but at the same time 



nl 



duj 



n=l 



(94) 



and therefore, using (83) and equating terms we obtain 



•••' w "^r ln ^.»J (95) 

as the general relationship between the normalized wavelet 
derivatives and the derivatives of the logarithm of the wavelet. 
Here B n (ci, C2, . . . c n ) is the nth-order complete Bell poly- 
nomial defined implicitly by (83). We then find 

*1;,8, 7 (W/J, 7 ) = (96) 

* 2;/ 3, 7 (w/3, 7 ) = -/3 7 (97) 

*3;/3, 7 (w/3 )7 ) = -07(7-3) (98) 

*4;/3, 7 (^, 7 ) = 3(0 7 ) 2 - 07 [(7 - 3) 2 + 2] (99) 

as first few values of the normalized wavelet derivatives at the 
peak frequency u)p tl . 

Appendix E 

Morse wavelet energy and peak frequencies 

For the generalized Morse wavelets, the wavelet "energy 
frequency" u5/? )7 defined in (12) is given by a ratio of gamma 
functions (43). Here we investigate why u3g i7 and u)p n should 
become indistinguishable for 7 = 3 and Pp n > 1, as was 
observed in Figure 4(d). 

Note that the asymptotic behavior of the gamma function is 
[24, eqn. 6.1.39] 

T (az + b) ~ V^e-^iazY^- 1 ' 2 , \z\ — > 00 (100) 

with |argz| < tt and a > 0; here f(z) <~ g{z), \z\ — > 00 



denotes lim 



1 



J(z)/g{z): 
r (a; + nr) 



a;(n-i)r r(x + r) 

for real and positive 2, n, and r. Choosing a; 
n = 2, one obtains 

1 r (20/7 + 2/7) 



1 as usual. It follows that 

- 1, |ar| — ► 00 (101) 
20/7 and 



W/J, 7 (20/ 7 ) 1/7 r (20/7 +I/7) 



1,0 — >oo (102) 



with fixed 7 but not, one may note, as 7 — > with fixed 0. 
Evaluating the left-hand side of (101) for n — 2 numerically 
(not shown), one finds that for r = 1/3, corresponding to 
the case 7 = 3, this ratio in fact remains very close to unity 
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for all x > 1 and rapidly approaches its asymptotic value as 
x increases. The minimum departure of the left-hand side of 
(101) from unity at a particular value of x is found to occur 
near r = 1/3 for all x > 1. Therefore the special properties 
of the 7 = 3 wavelets have their origins in the behavior of the 
gamma function ratio in (101) for r = 1/3. 
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